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ABSTRACT 


Prediction of demand and setting reorder levels by ex- 
ponential smoothing and mean absolute deviation is studied 
for an inventory system. The relation between the stand- 
ard deviation and mean absolute deviation of demand dis- 
tribution is discussed. Asymptotic first and second 
moments of the exponentially-smoothed forecast errors (MAD) 
are derived. Results obtained from simulation of several . 
normal systems indicated that the smoothing technique is 


inferior to the classical maximum likelihood estimation 


method. 


» SOHOOR 


33940 
TABLE OF CONTENTS 
Page 

Ti INTRODUCTION 9 
eT, RATIO OF A/o 14 
III. PROBABILITY DISTRIBUTION OF A 17 

AND OTHER ASSOCIATED STATISTICS 

A. EXPONENTIAL OF PAST OBSERVATIONS ey, 

(ESTIMATE OF THE EXPECTED VALUE OF 
DEMAND) 

B. FORECAST ERRORS 17 

C. JOINT DISTRIBUTION OF FORECAST ERRORS 19 

D. DISTRIBUTION oF |e, , | 23 
IV. MOMENTS or A 26 

A. MEAN VALUE OF A 26 

B. VARIANCE or A AG 
V. COMPARISON OF EXPONENTIAL SMOOTHING 39 

AND MLE BY COMPUTER SIMULATION 
APPENDIX A SIMULATION AT 
APPENDIX B NUMERICAL COMPUTATION OF VARIANCE oF A 50 
LIST OF REFERENCES 51 
INITIAL DISTRIBUTION LIST | 52 


DD 1473 = 











es 


IAA 


Ty. 


LIST OF TABLES 


Page 
VARIANCE OF A AND STANDARD 38 
DEVIATION OF 6 
SAMPLE CHARACTERISTICS 41 
PERCENTILES FOR REORDER LEVELS a3 
ACTUAL RISKS ATTAINED BY BOTH 46 


PROCEDURES FOR VARIOUS RISK LEVELS 





Qe 


TABLE OF SYMBOLS AND ABBREVIATIONS 


Demand during the ‘ia time period 
Exponential smoothing constant 

Forecasted demand during the ee period 
Correlation coefficient 

Risk level 

Mean absolute deviation of a probability 
distribution 

Mean demand 

Variance of the argument 

Expected value of the argument 

_th period forecast error 

Estimate of the standard deviation of the 
demand by maximum likelihood estimates 
Estimate of the standard deviation of the 
demand distribution by mean absolute 


deviation 





I. INTRODUCTION 


In recent years the technique of forecasting time 
series by exponential smoothing has received a considerable 
amount of attention in the literature. As a forecasting 
technique in inventory theory, it is particularly suited to 
the situation where demands are time dependent and corre- 
lated for much of its intuitive appeal rests on weighting 
the past in such a way that the most distant past receives 
less and less weight. When demands are not so related, 
smoothing may or may not be a suitable forecasting tech- 
nique. In particular, when demands are independent and 
from a common distribution such as for the so-called fast 
moving items, more classical methods like maximum 1likeli- 
hood may be preferable because of the optimal properties 
that are often demonstrable. Even in such circumstances, 
however, exponential smoothing is still a candidate for a 
forecasting scheme and it may very well be that, while not 
optimal , the difference in such criteria as stockout risk 
for example, may be negligible compared to optimal proce- 
dures from a practical point of view. In that case, fore- 
casting rules need not be changed for items differing in 
demand pattern and a certain uniformity may thus be 
achieved. 

The main purpose of this thesis is to investigate a 


model of an inventory system used by the Naval Supply 


System Command (NAV SUP) following Brown [Ref. 1]. The 
model, as described by Ref. 1 and Ref. 5, is briefly as 
follows: 

Single exponential smoothing is an operator given in 


Ref. 1l as 


t-1 
af k t 
s(x) =a ) BE x + Bex (1-1) 
k=0 
where x. is the total demand of the ie period 
® ais the smoothing constant with O.< @< 1] 


B=1-4 


Tntuitively speaking the smoothing operator is a 
method of forecasting by weighting past demands in such a 
way that the most recent demand is weighted by & and the 
remaining demands are given less and less weight as time 
progresses. 

The standard deviation, 0, of the demand distribution, 
which is needed to set the reorder levels, is estimated by 


the following three statistics: 


= eh ; 
e, = x, ~ S,_1 (x), t period forecast error, 
t-1 
A ak t 
A= a 
a ), 8 Fees + B be. (1-2) 
k=0 
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estimate of the mean absolute deviation of forecast errors, 


and 





o- fF Ata 
estimate of the standard deviation of the demand distribu- 
ron. 

This estimate of 0, as explained in Ref. 1, is based 
on the mean absolute deviation (MAD), A, of the forecast 
errors. Reorder levels, based on the predetermined risk 
level w (defined to be the probability of running out of 


stock), are found as follows: 
Ww = P(X > uw + ko) 


Here, Uw is the mean of X, and k is a constant which 
can be determined from the probability tables for a given 


risk level, vw. The reorder level is then given by 
Reorder Level = uw + ko (1-3) 


where k is the given above. 

In the analysis this model was slightly modified be- 
cause of the following considerations: 

dg Xoo the demand asked from the inventory system 
before it is placed in operation, is zero. Although one 
may predict a number for a fictitious inventory, there will 


be no observation for Xo° Hence ey 1S Zero. 


Hl 


as a’ and pN can be assumed to be approximately equal 
to zero for values of N > 30, 1£ 0.01 = @ = 0799. Ene 
of this, the modified form of (1-1) and (1-2) used in the 


analysis was 


t-1 
yak 
S.(x) =a ) BX, (1-4) 
k=0 
t-1 
B =a ) Ble, (ies 


k=0 


In this model, not only is demand forecasted by expo- 
nential smoothing but demand variability is estimated by 
means of exponentially smoothed absolute deviations. Such 
estimates are in turn used to set safety levels for inven- 
tory order rules. 

The model described above has been subjected to many 
criticisms. Zehna [Ref. 7] has indicated that the estimate 
of standard deviation might be very poor depending upon the 
underlying demand distribution. Based on the results ob- 
tained by Zehna [Ref. 8], the analysis in this paper takes 
on its starting point on investigation of the ratio A/o, 

a key element in the estimation of 0. It is quite evident 
that any weakness in the estimation procedure for 0 affects 
the reorder levels. An underestimate imposes higher risk 
and an overestimate causes unnecessary expenditures. In 
particular, the resulting consequences could be catastroph- 


ic ina military inventory system. 
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It is quite evident that the probability distributions 
of certain key statistics would be desirable in any exten- 
tion of the existing analysis. Although some of the dis- 
tributions required have been derived, the distribution of 
A, a most important statistic, could not be derived ex- 
plicitly. However, the first and second asymptotic moments 
Or A were obtained, and it is shown that the asymptotic 
variance of A is finite for the case where the demand has 
a normal distribution. Therefore, the next step of the 
analysis used simulation techniques to gain an insight into 
the behavior of the distribution of A, and to observe the 
effects of a finite variance. Reorder levels and their 
dispersions along with the risk levels attained are analyz- 
ed. The resulting errors due to the variability of safety 
levels is examined in detail and the consequences are 
evaluated from several points of view when demands are 
normal. Moreover the results are compared to maximum like- 
lihood methods establishing that a considerable difference 
exists in the two methods. For items satisfying the demand 
assumptions considered here, then, it is established quan- 
titatively that exponential smoothing is definitely inferior 


to more classical methods. 


aS 


DI g “RATIO OF AYa 


The mean absolute deviation of a probability distri- 


bution is given by 
A= s8([x - wl) . 


For a normal distribution the ratio A/o is the con- 
stant /2/7 = 0.8. Brown [Ref. 1] in his work on smoothing 
and prediction states, 

“Hence, the mean absolute deviation is proportional to 
the standard deviation, and the ratio may depend on the 
form of the distribution p(t), but dependence is slight." 

Now, as indicated by Zehna in [Ref. 7], this is not 
always the case. For example, consider a Poisson distri- 
bution with parameter ». A/o is then a function of A and 
it may be much smaller than 0.8 for small values of A 
[Ref.8]. However, there are some distributions where the 
approximation 0.8 for the ratio A/o is a good one under 
certain conditions. The quote above is a rather general 
statement which may be justified as follows: 

Let X be a continuous random variable with probability 
density function (p.d.f.) f(x) and let u and 1 be the upper 
and lower limits, respectively, of the random variable X on 
the real line. Then 


U 
d= e(x- ul) = [lx - ul tlxex 
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Or 


bs. 


U U U 
xEKX) Gx + i | f(x) dx + | xf(x)dx - yu | £ (x) dx. 
L L Ul U 


Hence, 
LL 
A= 2 | F(x) dx 
L 


Here, F(x) is the cumulative distribution of X,. 


Then 
u 
XY ee ae 
og oO pe tere: 
Let x = oy + ww. Then, 
m 
A 
m7 F(oy +) dy 
i= 
0 


Now suppose that F(x) is any distribution function such 


that 


F (x) eT 
oe 


Qle 
! 
8 


as 


where @(z) is the cumulative standard normal distribution 


function. Then for large values of u/o 


F(oy + wu) = Py) 


NES 


Hence 


QlIe> 
1R4 
NO 

Cy 

— 
O, 

N 
O, 
Ke 


Integrating by parts, the ratio then becomes A/o = J2/t. 
There are some distributions for which the above approxima- 
tion holds; for example erlang and chi-square. But for 
other distributions such a generalization may not apply. 


For instance, let X have the following mass function 


: if x= 
p(x) = 40.8 if x= 
Oy 1f x = 2a 
Then 
Nor=. O42a 
SRS Me 
and 


Ly Oe oe 


The error made by assuming A/o = 0.8 is 0.618, or 77% , 


Which is not negligible. 
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Itl. PROESEILISY DISTRIBUTION OF A 


AND OTHER ASSOCIATED STATISTICS 


Of central interest in the use of MAD is the ability 
to estimate the demand variability. The variance of the 
forecast error is always a function of the variance of de- 
mand, of. Therefore the MAD of the forecast errors is also 
a function of 0. Since MAD is a linear function of the 
standard deviation [Ref. 1] and if MAD of the forecast 
errors 1s estimated, then 0 can be estimated as described 
in [Ref. 1]. But in this estimation there are two points 
to be considered. 

i. Invariance Principle: 

Oo ine “are rellated to wach other by a aonstant, 
ue. , 

Ge = ca. 

But this relation is in the population, and it does not 
mean that the estimator of o and A will enjoy the same 
relationship necessarily. If these estimates mentioned 
above were MLE, then the invariance principle does hold. 
However, the definition of A indicates that this is not 
the case, 

2. The Constant of Proportionality of A/o: 

If the underlying demand distribution is not 
normal, then the assumption that the proportionality con- 


stant is approximately 0.8 needs to be validated. For this 


a / 


case the current model can be applied by a slight modifica- 


tion on 4G, that is, replacing /"/2 by the reciprocal of 


the appropriate proportionality constant. 


In this section we assume that demand is Normal (a,o*), 


where a is a constant. Derivation of the probability dis- 


an 
tribution of A requires probability distribution for the 
following statistics: 


ave EXPONENTIAL SMOOTHING OF PAST OBSERVATIONS 


(Estimate of the Expected Value of Demand) 


t-l 
SCX) =e ae Xk 
k=0 


Since all X's are mutually independent, the 
distribution of S (x) for large values of t is known to be 


a normal distribution with mean a and variance a7 (1 + Ta 
Bee FORECAST ERRORS 


eri has a normal distribution with the follow- 


ing mean and variance: 


uw = O for large values of t; 


= 2 


V ( = (0) hagas 


Sei 


for large values of t. 
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oe, JOINT DISTRIBUTION OF FORECAST ERRORS 
Let U and V be t-dimensional vectors defined 


as 


D 
(Tr 

Gq 
Gir 


an 


I 
~ 
I 
~ 


S 
il 


0) 

NO eeeoee0aee CT 
G Ga 

NJ eeeee0e0ee cr 


0) 
~ 
‘ i iaae a < 
f= 


where Uy. Be acy UL are constants. Let Y be defined as 


yY=uwu'vV 
or 
E 
Y = ), Uses 
_— 


Substituting e, in terms of the x's, 


c 1,=2 
= - = 
YX ), UK - L B Seen 
i=] k=0 
or 
E t 1-2 
“ieee Ux, ee ) ) ume x 
niieab a i-k-l 
oe 7k —O 


ES 


Hence, 


=! 


i-l 
a Rae alk ) 
sas Caan ici L B reaye 
k=0 : 


KK 
ll 
I Pet 


1=0 


Since the coefficients of the X's are constants and the X's 
are mutually independent normal variables, it follows that 
Y has a normal distribution. Hence by the definition of a 
multivariate normal distribution [Ref. 2], the forecast 
errors have a multivariate normal distribution. The corre- 
lation coefficient between ra and cae where 1 < j, can 
be computed as follows. 


Let 


NS aE ee where i1<j)j<t. 


on 


vee. = V(e._4) + Seis Ee), + 2Cov(e,_, e 


ue 


Cov(e,_.,e ) = + | v(y) = V(e,_3)-Vle,_.) | (3205) 


t-J 


Rewriting in terms of the X's, 


)-2 
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Rearranging terms on the right-hand side, 


jJ-1-2 
_ a ee a)-L-) v 
aa a + (1-a8 Xo; 
k=0 

aaa 

j-i \ wae 

Gare OR PR. 1, 

k=0 


Since the terms in the above expression are mutually inde- 


pendent, the variance of Y for large values of t can be 


found as follows: 


2(4-i-1) 4~i,2 
V(Y) ee of Ga + a2 ame a (a — op) ~ 1- i ai a* (+6) 
1-8 1 - B 


Swisetitutang V(Y), V(e,_3), and Vr.) into (3-1), 


a ca ee Saaee' 2 
Cove(e,_.-e, f aa eee T+h aS p24) * 1) (146) Ty) + 


eee ee 
aap?) - a) 


Hence, the asymptotic correlation coefficient between epi 
and e, . becomes 
t-J 
h JtE ( BO ey,2 (j-i-1) 4-i,2 
p=“ (1 + 7g (1-8 are =) rr 
j-i-1,2> 
(1-a8 yo )- 1 (3-2) 
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Now consider the following function, 


= if Gl + Seen el a-Y oe (1+8%) “) + (1-ag¥1)*>)_ 1 


Here y is a continuous variable. Then 


2y ia 
ay 7 & (= Aha log B+ 20% Wee ae OCP Loe, B)+ 


4 2 Bey , 
EE (- 22 log 8 + Ep 10g 8 ) 


Rearranging the terms, 


OL a apr 1+8)a@ | a (1+8) BY 
g2'= pY 1og a (- oe agg + UAL 
or 
a ie peeton e (ow 7 USE)o) = 7 apy log 8 
ay 28 


Now log & is negative; therefore, the first derivative of 
o', with respect to y is positive. Hence o', is a monotonic 
increasing function of y. Then positive integral values of 
y also lie on this same continuous curve. Actually, PE co- 
incides with 0', whenever y is a positive integer. Hence ¢ 
is a monotonic increasing function of j-i. The supremum of 


Oo is found by letting j-i ~" in (3-2): 


Sup (¢ Fd + +03 (1+1) + 1 )- 1 
0 


ZZ 





The minimum value of 0 is found by letting j-i = 1: 


eo. = Cat SQ ((1-1) + (148)7) + (1-097) - 1 


min 1+8 


= aa Ga + (14+B)a + BY ) = 3 


~it8 _j 
2 
The minimum value of 8 may be zero; hence, 0 = -0.5. 


min 


D. DISTRIBUTION OF jp 


Rewriting le | as 


t-1 


| | 7 
en "A Joe ae 


The term 


D 
hi 


2 oF 
o (1+ T+A) 


has, asymptotically, the chi-square distribution with one 
degree of freedom [Ref. 3]. 


The term 


Ze 


has, asymptotically, the chi distribution. Hence the mean 


value and variance of e 1S given by 


t-1 


E( le, 1) = ose /¢ (3-3) 
eee = sae (3-4) 


An attempt to derive the probability distribution of A was 
made utilizing the change-of-variable technique, and by 


defining the following transformation: 


t-1 
Bee) Si | 
= ere ee 
k=0 
Si ee ei 
ee eo, 


Ee ark 
5) am ce Yio 
: t-1 
Z _\ gk 
e. = + (A ), 8 (+ a) 
k=l 
This transformation is obviously not one-to-one; there 


are at sets in the domain of the transformation that are 
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transformed into a single set in the image. For all cases 


the absolute value of the Jacobian of transformation is 


one. The joint density of A, Vie anes Yeo then becomes 
he t-1 
f(y), : ae Y,_-17 9) = ? J. 2Vy" LY 5s oo. +(6-& i, B (4y,,) ) 
i=l k=1 


where g(x,, Kor cece x) is the p.d.f. of a multivariate 
normal distribution, and for every i one combination of the 
Signs is valid. The distribution of A then naturally re- 
quires the (t-1)-fold integration of the joint density of 
A p Vineees Ye - Although every parameter in the joint 
density is known, the labor required to arrive at an ex- 
plicitly-defined probability distribution makes the task 
very difficult. Therefore, the analysis concerning the 
dispersion of Oo was not completed because of the intracti- 


belity. 


ap 


A 
IV. MOMENTS OF A 


The next step of the analysis, then, was a simulation 
of the inventory system in order to gain an insight concern- 
ang the diecr Souci on Or A. It was seen that the reorder 
levels were fluctuating around the theoretical reorder lev- 
els, even after 1000 observations. This was related to the 
asymptotically finite variance of A as indicated by [Ref.7]. 
Although computer simulations implied variance of A does 
not vanish, the moments of A were derived in order that an 
analysis of the simulation results might be made. The mo- 
ments of A were derived for large values of t and are, 


therefore, asymptotic. 


A 
Ne MEAN VALUE OF A FOR LARGE VALUES OF t. 


t-l 

A oc 

A= @ eens abe = 
k=0 


ill 
© 


tol 
O 1s E ( le, | ) 
k=0 


(A) 


f 


Jad za ° 


“aS 
B. VARIANCE OF A 


, C-1 t-2 tol a 
v (A) = a7 ) py (| Soci +207 ) y) B* Joov (| 7 oe le = 
k=0 i=0 j=itl 
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covite,_;|. le,_,]) = Etle,_,e,-.1 - Bile, lee. .)) 


(4-2) 
Thus, substituting (3-3), (3-4), and (4-2) into (4-1), 


“A — 
woo es a _l, 


t-2 t-1 
20:2 ) ) atts E(le e | ) (4-3) 
t-i t-j 
i=@ j=i+) 


a” 
This last expression for V(A) requires the computation of 


E(le,_je,_5!)- The following derivations concern 


E(le,_je,_sl). 


t-1 


in Section III, is bivariate normal. The parameters of 


The joimt dustribution of e, . sand ep_4° 38 imcacated 


this distribution are 





Ou 

nal = E(e,_;) 0 o7 = 0 Js + I+A 
ai 

ae) - Ee) = 0 o.5 = 0 f1 + T+A 


and asymptotic © is given by (3-2). For notational con- 


venience, let 


Pe4| 


Then the joint p.d.f. 1S given [Ref. 4] by 


“4 
ee ne - ~™ <x < @ 


Z 
27 0, 0v1-p a eee 


where 
IL x xy y? 
oe Pe C FP oo i re 
QO ee y y 


In that case 


CO joe) 


e 
( | | xyl £(x,y)dy dx 


—_—O 


E(| xy} ) 


i 


= P es 
| |x| C-| y £(x,y)dy + | ¥ £ (x,y) dy ) dx 
- a O 


if) 


O O © 
| | xV © (tien Gyo. — | | xy £(x@idyawe 
7 22) —fe9) oO 


_— 


CO O oe) oe) 


| | my F(x, evox + | | xy £ Gy) chads 
oO -° Oo Oo 


However, by the definition of correlation coefficient, 


E (XY) 


O O O fore) 
| | xy £(x,y)dydx + { i Saye (oc paneny Cx 
—lee) —(e2) —les, O 


oe) CO 


co O | rm 
( | xy £(x,y)dydx + | | may Elo) Cy cx 
O -0 O O 
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Thus, 
O C 6 


O 
E(| xy| )+E(xXy) = 2 C | | xy £(x,y)dydx + | | xy £(x,y) dydx ) 


Making a change of variable in the first integral, let 


Z = -x 
u= -y 
then 
O O fo) foe) 
| | xy £(x,y) dydx = | zu £(-z,-u) dudz 
“7 "-2 oepaere: 


Since £(x,y) = £(-x,-y), it follows that 


O O foe) foe) 
| | xyf(x,y)dydx = [ | Zu £(z, ~) adindz 
a ome 
Thus 
e(}xv|) + E(xy) = 4 | | xy £(x,y)dyax (4-4) 
» © 


The right-hand side of (4-4), without the coefficient 4, is 





ea _ tyab)- 
Oo fo oe) 20° © 20° (1-p*) 
xe r ye 
| | xy E(x, yyeydx = | a — dy ax 
Oo oO O Coan I J20 o/1-p? 


(2=5) 
where b is the conditional mean of Y given X, and is given 
by [Ref. a, 


Le 
je. = 
Ge OX 
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The conditional standard deviation of Y given X is given by 








, 2 
Cua f Le 
y p 
Let 
So FE Hen 
G 
a (y-b) 7 2 
5 =D 5 a? 2° 
oo 20 co Z2¢ co 


x 
xe | ye See Ge or az 




















(4-6) 
Integrating the inner term first: 
2 ne 
foe) _ 42 5 ae, 
S292 2 7 Ge aa So aS (1 - &- 25 
“p fan J2T 
c 
where ®(x) is the cumulative normal distribution. Sub- 
stituting this expression into (4-6), 
2 2 P. 
x By * oe 5 
, 20.7 : 5 ; 20 7 . 2 
oe = F e 2e reas e = ; 
fate fe ae [x __ 
O oJ 2m oc 2A Oo ox 27 
2 a 
xX —_ 
i og 2 b 20 2 
ase e eb O- aa ‘ 
| ax 2 ax (4-7) 
O oO /2m O oj 2m 
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Denote the terms on the right-hand side as l, 2, and 3, 
respectively, and compute each one separately as follows: 


Ll. Suibstituting values of c and b into 1, 














2 
x 
2 2 ee a 
=e int 2 2 20,.*(1-0) 
> = 862 ep HG, og /1l-p*x e 
Ce a Me ae Le Se 
| 276 dx = 270 ox 
O x ee x 
ae 
o /1-0 
ae _o* 
7 27 a. (1 po”) 
9 3 
; O a, (1-9 )2 
- 27 
2. Substituting the value of b into 2, 
- x? x 
2 a 
20 2 
I ee a o oO 7° 5 a 
ee a 
Oo Os J2T wz PO «® 
Co OC QO 
_ aes 
2 
3. Substituting the values of c and b into 3, 
Z 
- gal 
- 20.7 P 
O si J/1-p2 y o Jor 05. 
x x 
x2 
ome ®) po - 2G: Ere ae et 
y 2; 
oy 326 (O 2 
x 
where p= ge oo 


Let. x/o,, = y in (4-8); then (4-8) becomes 

















2 
Oo /p y) =f /2 2 
ze | ye ( = dt dy 
J27 O —-~% /27 
where r= Py 
v1-( 
Integrating by parts, 
2 a 
oo -© = 2 
Z - 
du =-ye a way Va "yy. | =< aa ee 
=o S29 
2 a r* 
- + 7. 2 geen 
oa 2 =, = e p 
Le lee dv = adtmy dy 
j -~ ,/27 /1-p7 
Eq. (4-8) becomes 
2 
. 2 .-_ 
Yo ge ee 2 
i eG | Pye 2/1-p 
ae | eu Kk) OS ea «Cay + 
2 2 
ee Y_. \ eee ee 
IX e 2 e = 
‘ Ct scivyrae 
° Toe 
where 
Co 0 / 
1 oe 
J/27 
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The first term in the above equation is zero; hence, 


(4-8) becomes 





2 Es 

iL = is 

Pred p*(1-p*) 7 aay aa : Te 
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(4-10) 
This result, due to the fact that y changes over positive 
numbers and ~P is between -.5 and zero, makes the argument 


of ® (-) positive. Hence (4-10) becomes 


Z i 
teow. es dealing 
T 
5 i 5 + |e | e dt dy = M 
Let 
e— V J2 
Y = (8G 
then 
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Since the series is uniformly convergent, it follows that 
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It was previously shown that € is always between -.5 and 0. 


Therefore, 


ee ee | 
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Hence the series in the above expression is the expansion 


of Arctan C=? liseses Gl). - 
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Thus Eq. (4-10) becomes 





1 1 = 
5 ‘i 5 + Wen Arctan G ¥ ) (4-11) 


Substituting (4-11) into (4-9) yields 





1, 
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Combining the results for 1, 2, and 3 into (4-7), 











and sub- 
Stituting into (4-5) yields 
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But o.. and c., are both equal to @ i 1 see) - Hence, 
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)) (4-13) 
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Therefore, 


t-2 t-l 


v (A) ion 2 7 a a: o74 207 ) Pele seal) 


(1487) -(148)* i oe 


(4-14) 


where E(le | ) is given by (4-13). 


t-j Sti 


It is clear that to arrive at an expression for v(A) 
in the form of (4-14) many tedious and difficult calcula- 
tions must be made. In fact, in order to arrive at an 
explicit value for v(A), one must deal with an arctangent 
term. However, because o is common to each term, the 
ratio v (A) /o? can be defined in terms of the known con- 
stants, and an approximate numerical answer can be found 
for the ratio v(A) /o*. The computer program in Appendix 
B computes the first two terms in (4-13) exactly, and the 
third term approximately. Briefly, the computation of the 
third term was accomplished as follows: 

Eq. (4-13) is a bounded number for a given @ and 8. 
It cannot be made greater than some finite number d, what- 
ever be the value of p . (Note that p is the only varuamme 


in (4-13) that depends on i and j, and it is bounded.) 
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Therefore, in the double summation of (4-14), for some 
finite value of i and j, higher order terms will get smaller 
than a specified number m which, in this case, was chosen 


as von +? The program was written so that the higher terms 


10 were not included in the result, and 


smaller than 10 
were subsequently dropped from the summation. The upper 
bound of the error introduced by omitting the higher terms 


from the summation can be found as follows: 


E(| xy) is bounded from above by 


Since the summation is infinite, 


Upper Bound of Error 20 7. ) ai 
1=0 j=itl 


Aeerion.” B 


(1587) - 


a4 geile 


(1+A) 


First, this error is additive and hence is added to 
the result. Second, whatever be the value of 8, it is 
always much smaller than one. Hence, the results contained 
in Table 1, obtained by the computer program for various 


values of @ were almost exact. 
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SMOOTHING 
CONSTANT 


Table l. 


0.0204 
On0 S25 
0.0461 
0.0614 
0.0785 
OREO See, 
Ore 
0.1440 
Os lo 
0. 2026 
0.4267 


“A 
Variance of 4 
Oo 


VARIANCE OF (A) /o? 


and Standard Deviation of 
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STANDARD 
DEVIATION OF 


0.1745 
©. 2a7s 
gh a) 
@.2905 
O25 2 
O. 3568 
0.3875 
0.4187 
0.4495 
0.4802 
0.6341 
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V. COMPARISON OF EXPONENTIAL SMOOTHING 


AND MLE BY COMPUTER SIMULATION 


For the reason explained previously, simulation was 
utilized to explore the parts of the model where theoreti- 
cal work was very difficult, and was also used to compare 
the exponential smoothing with MLE. The reason that MLE 
was chosen as the alternative procedure was because it has 
the property of an asymptotically vanishing variance. As 
indicated in [Ref. 8], the maximum likelihood estimator, 


0 , of O is given by the formula 


where x is the sample mean. The statistic @ has a chi- 


square distribution [Ref. 3]; ain fact, 
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From these formulas, 
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The computer program used for the simulation is included in 
Appendix A. The simulation procedure and ensuing results 
are as follows: 

One thousand normally-distributed random numbers, re- 
presenting the demand, were generated by a pseudorandom 
number generator. Various values of the parameters, JU, OG, 
and w were used. For every triplet of the parameters, re- 
order levels, risk levels attained, estimates of O and 
sample mean, and sample standard deviation of G and of the 
estimate of o0 by MAD were computed. The value of & was 0.2 
for all cases. Table 2 presents the results achieved by 
the simulation. As Table 2 indicates, the theoretical 
values of the standard deviations of both estimates of 0 
were in complete agreement with the sample standard devia- 
tions of these same estimates. The variance of @ is higher 
than that of @ in every case. The last two columns ae 
the mean square error (MSE) of O and G. The MSE is given 
by the variance plus the square of the bias. Although ao 
is biased for 1000 observations, its bias can be neglected. 
Then the MSE of 0 and G are determined by their variances. 
The comparison of the MSE columns are quite striking; the 
MSE for 6 reaches values as high as 100 times that for . 
By noting the similarity between the theoretical values and 
the simulation results, the effect on reorder levels was 
analyzed as follows: 

Computation of reorder levels was made by utilizing 


Eq. (3-3). First, the theoretical values of the reorder 
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levels were computed for every fixed parameter triplet u, 
o, » . Then the reorder levels, found by smoothing and 
MLE, were compared with the theoretical values. This com- 
parison is summarized in Table 3. The value Q) in Table 4 
is the percentage (rounded off to the nearest integer) of 
1000 reorder levels found to be within one unit of the 
theoretical reorder level. The first column under each of 
Pros and Po. lists the results for smooth- 


ed estimates; the second column under each heading lists 


the headings Qu. 


MLE. The quantity P represents the number of units about 


50 
the theoretical reorder level within which 50 per cent of 
the computed reorder levels were found; Pos is a similar 
computation for 95 per cent of the reorder levels. Table 
3 clearly illustrates the superiority of MLE compared to 
exponential smoothing. For the case where smoothing was 
used to estimate 0, reorder levels fluctuated around fhe 
theoretical reorder level, causing either an excessive 
purchase, or imposing more risk than desired. The magni- 
tudes of these fluctuations are clearly evident from the 
table 

Table 4 indicates the risk levels attained by both 
methods at the 1000th observation. For every predetermined 
parameter pair w and o, five different risk levels were 
selected. For every triplet, the theoretical reorder 
levels were compared with the computed reorder levels of 


both procedures. The actual risk levels attained were then 


computed. 
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PARAMETER 
TRIPLETS 


40,4,.01 
40,4,.11 
40,4, .50 


5055, .01 
a0, , «lel 
50), 5,.'50 


LeOe tO, 62 
0710, a2 
awe, tO,.50 


ZOO 20ije20 1. 
200,20,-L1 
200 202 50 


200,30, .0l1 
a0, 30% . 11 
S007 30,.. 50 


400,40, .01 
400,40, .11 
400,40, .50 


moo, 50, .01 
506,50, .11 
ow, 00, .50 


600,60, .01 
S00760,.11 
600,60, .50 


7OO,70,.01 
woo, /0,.L1 
730,70, «50 


600,50, .01 
800,80,.11 
800,80, .50 





Table 3. Percentiles for Reorder 
Levels. 
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PARAMETER 
TRIPLETS 


9008 907.02 
90GB 07. 2 Sia 29 





1000,100,.01 }| 1 25145 3 
LOCO. OO. i 1 8 | 29 
LeGO; LOGE. 50 1 6 4; 25 2 
Table 3.(Continued) ~- Percentiles for Reorder 
Levels. 
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As Table 4 indicates, the risk levels attained by MLE 
are found to be very close to the desired risk levels. 
However, the risk levels of the smoothing method were 
found to fluctuate over a wide range. 

These results, gained by the simulation, indicated 
that the MLE technique is to be preferred over the tech- 
nique of exponential smoothing. However, due to the lack 
of theoretical analysis concerning the distribution of 
these results are not considered conclusive. However, this 
analysis does indicate that whenever a constant mean de- 
mand can be justified, the MLE technique is to be prefer- 


red to the exponential smoothing technique. 
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